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Abstract 

The goal of this project was to investigate the effects and possible developmental disease implication of chronic dietary 
TCDD exposure on global gene expression anchored to histopathologic analysis in juvenile zebrafish by functional genomic, 
histopathologic and analytic chemistry methods. Specifically, juvenile zebrafish were fed Biodiet starter with TCDD added at 
0, 0.1, 1, 10 and 100 ppb, and fish were sampled following 0, 7, 14, 28 and 42 d after initiation of the exposure. TCDD 
accumulated in a dose- and time-dependent manner and 100 ppb TCDD caused TCDD accumulation in female (15.49 ppb) 
and male (18.04 ppb) fish at 28 d post exposure. Dietary TCDD caused multiple lesions in liver, kidney, intestine and ovary 
of zebrafish and functional dysregulation such as depletion of glycogen in liver, retrobulbar edema, degeneration of nasal 
neurosensory epithelium, underdevelopment of intestine, and diminution in the fraction of ovarian follicles containing 
vitellogenic oocytes. Importantly, lesions in nasal epithelium and evidence of endocrine disruption based on alternatively 
spliced vasa transcripts are two novel and significant results of this study. Microarray gene expression analysis comparing 
vehicle control to dietary TCDD revealed dysregulated genes involved in pathways associated with cardiac necrosis/cell 
death, cardiac fibrosis, renal necrosis/cell death and liver necrosis/cell death. These baseline toxicological effects provide 
evidence for the potential mechanisms of developmental dysfunctions induced by TCDD and vasa as a biomarker for 
ovarian developmental disruption. 
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Introduction 

TCDD is a highly toxic, widespread environmental chemical 
that is a by-product of incinerating plastics, chlorinated industrial 
chemicals and other hydrocarbons. TCDD has been shown to 
persist in the environment, bioaccumulate in humans and wildlife 
through trophic transfer [1] and contributes to adverse disease- 
specific effects such as cancer and cardiovascular toxicity [2-5]. 
Zebrafish, as an aquatic animal model, continues to play a central 
role in TCDD toxicological developmental and genetic research 
[6]. TCDD has been found to cause precursor and disease- 
associated adverse effects in zebrafish embryos, such as craniofa- 
cial malformations in jaw, edema in yolk sac and pericardium, 
decreases in cardiac output [7], spinal deformity [8], disruption of 
development of common cardinal vein [9] , and a developmental 



reduction in blood flow in mesencephalic vein [10]. In TCDD- 
treated adult zebrafish, collections of fully developed disease 
phenotypes and dysfunction have been documented. For example, 
TCDD has been found to inhibit caudal fm regeneration in 
zebrafish, suggesting dysfunction of tissue's self-repairing abilities 
by its inhibition of cellular differentiation and proliferation [1 1]. In 
the reproductive system, TCDD represses the follicle maturation 
in ovarian development and inhibits estradiol biosynthesis [12]. 

In spite of these wide-ranging and impactful early embryonic 
dysregulation findings and fully developed adult disease implica- 
tions, there are few studies examining the onset and development 
of the dysregulation in fish during the critical developmental 
transition from embryo to adult. For example, neither toxicolog- 
ical nor toxicogenomic effects of dietary TCDD exposure on 
juvenile zebrafish are well documented. Although it is well 
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understood that key dysregulation critical to adult health are 
initiated in the embryonic stage, many diseases and exposure- 
pathology emerge during the juvenile stages (e.g. so called "early" 
onset diseases) or precursors to the adult problems are not 
observed in embryonic studies may be identified in the juvenile 
stage. There is a gap in our understanding of the pathology, gene 
expression and dysregulation of juvenile zebrafish due to dietary 
TCDD exposure. We hypothesize that early onset precursors to 
adult diseases and toxicogenomic dysregulation can be detected in 
juvenile fish and that juveniles may exhibit dysregulation not 
observed in either embryos or adults. 

The objective of this study is to create a baseline of evidence of 
the phenotypes and genetic dysregulation leading to a better 
understanding of the molecular mechanism (s) involved in the 
onset of TCDD-induced disease exhibited in zebrafish. Chronic 
dietary TCDD exposure was conducted by exposing juvenile 
zebrafish to dietary TCDD at 0, 0. 1, 1, 10 or 100 ppb TCDD for 
42 d. Accumulations of TCDD in fish were measured in whole fish 
bodies. TCDD-induced toxicities were investigated using histo- 
pathologic analysis, and global gene expression in fish treated with 
TCDD was analyzed through microarray experiments. TCDD 
accumulation was dose- and time-dependent in zebrafish, 
although no significant mortality was observed even at accumu- 
lated doses up to 10 ppb. Juvenile zebrafish showed pathologic 
alterations due to TCDD exposure at 28 and 42 d, especially in 
disruption of ovarian development and degeneration of nasal 
epithelium. Functional analysis of microarray data from TCDD- 
treated fish revealed a correlation between gene expression and 
histopathologic outcomes. 

Methods 

Ethics statement 

Zebrafish were handled following protocols approved by the 
Animal Care and Use Committee (ACUC) of the University of 
Wisconsin-Milwaukee (UWM). 

Chemicals and diet 

TCDD was purchased from the National Cancer Institute- 
Chemical Carcinogen Reference Standard Repository, Midwest 
Research Institute (Kansas City, MO) and stored at room 
temperature. Dose formulations for treatment were prepared by 
dissolving TCDD in ethanol and mixing the appropriate amount 
into Biodiet starter #2 (Bio-Oregon, Inc., Warrenton, OR), using 
1 L ethanol per kilogram of diet. The ethanol in the mixture was 
evaporated under a fume hood to remove the vehicle. 

Animals 

Juvenile wild-type zebrafish (stock originally purchased from 
EkkWill Waterlife Resources, Ruskin, FL) (0.16 g, 2.10 cm) were 
reared in 28°C dechlorinated, filtered municipal water (DFMW) in 
the NIEHS Children's Environmental Health Sciences Core 
Center (CEHSC) animal facility. Fish from a single treatment were 
housed in a single flow-through 40 L polycarbonate tank to reduce 
hazardous waste. Zebrafish were fed Biodiet starter (4% body 
weight per day) for 42 d with TCDD added at 0 ppb, 0.1 ppb, 
1 ppb, 10 ppb or 100 ppb (nominal concentration, ethanol used as 
vehicle). Each treatment had one tank with 300 fish (N = 300). 
Throughout any animal treatment, fish were visually monitored 
three to four times daily as part of normal husbandry. Any 
individual that exhibited signs of distress (the inability to maintain 
an upright position, holding fins close to body for a prolonged 
period, erratic swimming, remaining close to the bottom or top of 
the tank, inflammation around the gills, excessive and/ or difficult 



respiration, or suppression of pigmentation) was immediately 
humanely euthanized. For euthanasia, fish were anaesthetized in 
100 mg/1 of buffered MS222 (tricaine mesylate) solution, then, 
killed by cervical transection or immersion in liquid nitrogen. Fifty 
fish were collected at 7, 14, 28 or 42 days (d) respectively, following 
initiation of TCDD diet, and 15 of them were measured for body 
weight and length. Samples were stored at — 80°C until used for 
experiments. 

TCDD assessment 

All samples have been analyzed by the French National 
Reference Laboratory for PCDD/Fs and PCBs in food and feed, 
using the method described below that was accredited according 
to the ISO 17025 standard. Details of the sample preparation 
method can be found elsewhere [13-15]. 

Five individuals of each group were pooled for extraction using 
a Pressurized Liquid Extraction system (ASE, Dionex, Sunnyvale, 
CA, USA). Three successive static extraction cycles (5 min each) 
were performed using a mixture of toluene/ acetone (70:30, v/v) at 
100 bar and 120°C. Extracts were aliquoted before purification 
steps, according to expected TCDD content. Clean-up and 
fractionation of TCDD was carried out using liquid chromatog- 
raphy columns with sulphuric acid silica, Florisil and carbon 
(Carbopack C) in that order. 

AU samples were analyzed by gas chromatography coupled to 
high resolution mass spectrometry (GC-HRMS). Identification 
and quantification of TCDD was achieved using the isotope- 
dilution method. ^^Ci2-labeled TCDD was used as internal 
standard. GC separation was performed on a DB-5MS capillary 
column (30 mxO.25 mm, 0.25 |im) in HP 7890 instrument (HP, 
Palo Alto, CA, USA). HRMS measurements were achieved on a 
JMS 800D electromagnetic instruments (Jeol, Tokyo, Japan), 
operating at a resolution of 10,000 in the selected ion-monitoring 
acquisition mode after Electroionisation (42 eV). 

Histopathologic analysis 

Five to ten fish were collected after 28 and 42 d of dietary 
exposure from each treatment group for histopathologic analysis. 
Briefly, fish were anaesthetized in 100 mg/1 of MS222, and then 
sacrificed by cervical transection. The coelomic cavity was opened 
allowing perfusion of internal organs. A transverse incision at the 
level of the anterior spinal cord enhanced fixation of brain and 
kidney. Tails were removed, and carcass was placed into plastic 
cassettes separately. The samples were fixed in Davison's solution 
(150 ml ethanol, 95%; 100 ml formaldehyde 37%, 50 ml glacial 
acetic acid and 58 ml distilled water) for 24 h, rinsed in cold tap 
water (24 h), decalcified in Cal-EX for 24 h (Fisher Scientific, 
Pittsburgh, PA), and rinsed again in cold tap water (24 h). 
Dehydration occurred in 50% and then 70% ethanol for 24 h 
each. The samples were subsequently embedded in paraffin. 
Sagittal-step sections were cut near midline and medial to the eye 
(2-4 per fish), mounted on glass slides and stained with 
hematoxylin and eosin by Mass Histology Service, Inc. (W orcester, 
MA). Histopathologic analysis was performed at Oregon State 
University (Corvallis, OR) by a veterinary pathologist certified by 
the American College of Veterinary Pathologists, having 30 years 
of experience in fish pathology research and diagnostic pathology. 

RNA isolation 

Total RNA from individual fish was isolated using Trizol 
reagent (Invitrogen, Carlsbad, CA). One ml of Trizol was used per 
individual fish, which were then homogenized with Power Gen 
500 (Fisher Scientific, PA). Chloroform (0.2 ml) was added for 
dissociation of nucleoprotein and mixed by shaking the tubes by 
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hand for 15 s. The samples were then centrifuged at 12000 g for 
15 min at 8°C. The upper aqueous phase was moved to a clean 
1.5 ml microcentrifuge tube, and 0.5 ml of isopropyl alcohol was 
added for precipitation of the RNA by incubating the samples at 
room temperature for 10 min. The samples were centrifuged at 
12000 g for 15 min at 8°C. Then the supernatant was removed 
and 1 ml of 75% ethanol was added to wash the pellet, and the 
samples were centrifuged at 12000 g for 5 min at 8°C. The 
samples were air-dried and then dissolved in RNase-free water. 
The concentration of RNA was determined with a ND-1000 
NanoDrop Spectrophotometer (Thermo Scientific, DE). 

RNA purification was performed using the RNeasy MiniElute 
cleanup kit (Qiagen, CA) following the manufacturer's protocol. 
Pre-cleaned total RNA samples (<45 |ig in a total volume of 
100 |il RNase-free water) were digested with 2.5 |il DNase I 
(1500 U) for 10 min at room temperature. Digested RNA was 
eluted through a RNeasy MinElute Spin Column with 350 |il 
RET buffer and 250 |il 100% ethanol to bind RNA to the 
membrane of the column. The column was washed with 500 |il 
80% ethanol after being washed with 500 |il RPE buffer, and 
RNA was eluted in 20 |il RNase-free water. The concentration of 
RNA was determined using the ND-1000 NanoDrop Spectro- 
photometer. Quality of zebrafish total RNA was determined with 
an Experion RNA StdSens chip (Bio-Rad, CA). 

Sex designation 

Zebrafish at the juvenile stage were placed into the TCDD 
feeding experiment and typically achieve gonadal maturation after 
3-4 weeks. The body cavity of each fish was opened and the 
presumptive sex determined by visual assessment of the gonad. 
Unexpectedly, all fish in the 100 ppb-treated group were lacking 
mature or undifferentiated gonads at 42 d. 

We felt the need to determine the sex of each individual fish to 
account for differences in response between males and females. 
Since visual sex-determination is not 100% reliable, a novel RT- 
PCR-based sexing method was developed and validated. It has 
been reported that male and female zebrafish have two different 
splicing variants of the gene vasa after 25 dpf that differ in length 
depending on inclusion or exclusion of exon 4 [16]. Male zebrafish 
mainly express short vasa mRNA without exon 4 and have very 
low expression of long vasa containing exon 4, and female 
zebrafish have both short and long vasa mRNA splice variants 
(Figure 1). In order to distinguish male and female by RT-PCR 
based on the numbers of bands of RT-PCR products on the gel, a 
pair of primers was designed across from exon 2 to exon 6 as 
shown in Figure 1 . Both visual observation of the gonad and PCR 
with vasa primers were used to sex TCDD-treated zebrafish 
(Table 1). 



Quantification of exon-4 transcripts of vasa mRNA 

In order to identify presumptive TCDD-induced disruption of 
splicing of vasa mRNA between male and female zebrafish, 
another pair of primers was designed for testing levels of exon-4 
transcripts in two qRT-PCR experiments: 1) sexed male and 
female zebrafish in all treatment groups at 28 d; and 2) sexed 
male, female and immature zebrafish in 0 and 100 ppb groups at 
each collection time point. 

For quantification of the exon-4 transcripts in qRT-PCR 
experiments, RNA from female and male fish in 0 ppb groups at 
28 d was pooled and used as a common reference for the 
comparison. We identified guanine nucleotide-binding protein (G 
protein), beta polypeptide 2 -like 1 {gnb2ll) as a normalizer gene 
based on this and previous studies of TCDD-treated zebrafish in 
this laboratory. The sequences of primers for gnb2ll are: 
gnb21 1 _forward: 5 ' -AGCTGAGGC AGGAC ATC ATT-3 ' ; 

gnb21 l_reverse: 5 ' -GCTTTATCTGGTTCCGATGG-3 ' . 

Regarding QPCR, first-strand cDNA was synthesized from 
1 |ig of total RNA from individual fish with oligo-dT using an 
Affinity S crip tTM Multi Temperature cDNA synthesis Kit (Agilent 
Technologies, Inc. Santa Clara, CA) following the manufacturer's 
protocol. PCR amplification was performed with the Mx3000P 
QPCR Systems (Agilent) in 10 |il reactions using 1 |il of cDNA 
(10 ng of total input RNA), 200 nM each of forward and reverse 
primer and 1 x Power SYBR Green PCR Master Mix (Applied 
Biosystems, Carlsbad, California). The real-time PCR program 
consisted of 1 cycle of 95°C for 9 min; and 40 cycles of 95°C for 
15 s, 59°C 30 s and 72°C for 30 s, with fluorescence detection at 
the end of each 59°C step. The products were determined by 
melting curves. Relative gene expression data (fold-changes) were 
calculated using the mathematical model described by Pfaffl [17]. 
Statistical analysis of QPCR data was performed using SigmaStat 
11.0 (Systat Software, Inc., CA). One-way analysis of variance 
(ANOVA) was used to detect the significances of treatments 
compared to control group at 28 d. 

Microarray experiments 

We used individual fish to perform the microarray experiments 
and assess biologic variability in gene expression. Two microarray 
experiments were designed based on time, dose, and designated 
sex (Figure 2): (1) Dose-dependent experiment: gene expression 
testing conducted in TCDD-treated (at 0, 0.1, 1, 10 or 100 ppb) 
for male and female sexed zebrafish at 28 d; and (2) Time- 
dependent experiment: gene expression testing of 100 ppb-TCDD 
treated male and female sexed zebrafish at 7, 14, 28 and 42 d. 

NimbleGen Gene Expression 12x135 K zebrafish microarrays 
were used for genome-wide expression analysis of TCDD-treated 
zebrafish. RNA samples from six zebrafish (three male and three 
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Figure 1. Two vasa splicing variants in zebrafish (Krevel eta/., 2004). Male zebrafish have short vasa mRNA and female zebrafish have both 
short and long vasa mRNA. A pair of primers was designed crossing exon 2 and exon 6 (arrow). The sequences of primers are: vasa_long_forward 5'- 
CCCAATATGGATGACTGGGAG-3'; vasa_long_reverse 5'-CAGCACCTCCTGTATAAAAGC-3'. A PCR method was applied to distinguish male and female 
based on the numbers of bands of PCR products on the gel. Female fish have two bands and male fish have one band. In addition, expression of exon 
4 within vasa gene was quantified by real-time QPCR. The sequences are: exon4_forward: 5'-CTTCACGCGTGTCCACCTG-3'; exon4_reverse: 5'- 
GCTCTTG A ATCCTCC ATC AG A ACC-3 ' . 
doi:1 0.1 371 /journal.pone.01 0091 0.gOOl 
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Table 1. Sex designation of dietary TCDD-treated zebrafish using both observation and QPCR method. 





Visual observation 


PGR by vasa 


Designated sex 


testis 


male 


male 


ovary 


female 


female 


Abnormal gonad 


female 


female (immature) 



doi:1 0.1 371/journal.pone.01 0091 0.tOOl 



female, when available) were randomly selected. cDNA was 
synthesized using the Superscript Double-Stranded cDNA syn- 
thesis Kit (Invitrogen, Carlsbad, CA). Ten |ig total RNA from 
individual fish was denatured at 70°C for 5 min with oligo-dT 
primers, and then was reverse transcribed to the first stand cDNA 
at 42°C for 60 min with SuperScript II enzyme. The second- 
strand cDNA was synthesized at 16°C with DNA ligase, DNA 
polymerase I and RNaseH for two hours. After digestion with 
RNase A solution, the double-stranded DNA was cleaned in a 
Phase Lock tube (Fisher Scientific, PA) with phenol: chloroform: 
isoamyl alcohol (25:24:1, Ambion, NY), and then precipitated in 
ice-cold ethanol with 7.5 M ammonium acetate and 5 mg/ml 
glycogen, and dried to a pellet in a DNA 1 20 SpeedVac (Thermo 
Scientific, DE). The DNA was then cleaned using QIAquick PGR 
purification kit (Qiagen, CA). The quality of DNA was determined 
with an Experion RNA HighSens chip (Bio-Rad, CA). One |ig 
double-stranded DNA was used for Cy 3 -labeling reaction with 
random primer and Klenow fragment (3'->5' exo-) using a One- 
Color DNA labeling Kit (NimbleGen, WI). Cy3-labeled DNA was 
precipitated in ice-cold ethanol with 5 M NaCl and dried to a 
pellet in a SpeedVac, then rehydrated in nuclease-free water. The 
concentration of Cy3-labeled DNA was determined using a ND- 
1000 NanoDrop Spectrophotometer. Cy3-labeled double stranded 
DNA was hybridized to a NimbleGen 12x135 K array using a 
NimbleGen hybridization kit; the hybridization was performed in 
the NimbleGen Hybridization System (NimbleGen, WI) at 42°C 
for 16 h. Arrays were washed with NimbleGen Wash buffers and 
1 M DTT, and then spin dried in a NimbleGen microarray drier. 

Hybridized microarrays were sent to MOgene, EC (St, Eouis, 
MO) for scanning. The arrays were scanned at 2 |im resolution on 





0.1 


1 


10 


100 


0 


7d 




*3/3/3 


3/3 


14 d 


*3/3/3 


3/3 


28 d 


3/3 


3/3 


3/3 


3/3 


3/3 


42 d 




*3/0/3 


3/3 



Dose-dependent 



I 

Time-dependent 



Figure 2. Microarray experiment design for TCDD-treated 
zebrafish. Zebrafish were fed with TCDD at 0, 0.1, 1, 10 and 
100 ppb TCDD in food for 42 d. Fish were collected at 7, 14, 28 and 
42 d. Two microarray experiments were conducted. Dose-dependent 
experiment: fish from each treatment at 28 d for microarray analysis; 
Time-dependent experiment: fish from 0 and 100 ppb groups at 7, 14, 
28 and 42 d for microarray experiment. Three female and three male 
fish from each group were used for microarray experiments with 
NimbleGen 12x135 K zebrafish microarray system was used for this 
project. The asterisk (*) means three immature fish in 1 00 ppb groups at 
7, 14 and 42 d were used for microarrays, and no mature female fish 
were found in 100 ppb group at 42 d. 
doi:1 0.1 371/journal.pone.01 0091 0.g002 



a NimbleGen MS200 scanner with auto-gain adjust. The TIFF 
images were gridded and signals extracted using NimbleScan v. 
2.6. Expression data were normalized with the Robust Multichip 
Average (RMA) algorithm as described previously [18]. 

Microarray data analysis 

Differential expression analysis was performed with Multiple 
Array Viewer (MeV, Dana-Farber Cancer Institute, MA) software, 
version 4.8.1 [19,20]. Three statistical methods were reviewed and 
tested for analysis: Bayesian Estimation of Temporal Regulation 
(BETR), one-way ANOVA, and Significance Analysis of Micro- 
arrays (SAM). False discovery rate (FDR) was controlled with the 
Benjamini and Hochberg procedure [21]. For both experiments 
and all expression data, BETR consistently predicted the lowest 
FDR of the three methods, and was therefore selected for 
differential expression analysis. Statistical significance was set at 
FDR <0.05. We applied K-means clustering using MeV software 
to cluster genes based on their expression profiles across doses and 
time. 

Gene Ontology (GO) annotations (categories of biological 
process [BP]) were achieved by performing BLASTX between 
array sequences against all protein sequences of the Amigo 
database at amigo.geneontology.org [22]. Zebrafish sequences 
were extracted to align against the Amigo database based on their 
accession numbers on the array from different sources, including 
the NCBI (http://www.ncbi.nlm.nih.gov/UniGene/), VEGA 
(http://vega.sanger.ac.uk/Danio_rerio/) and Ensembl (http:// 
www.ensembl.org/Danio_rerio). The "GOHyperGAll" script 
from BioConductor was employed to conduct GO enrichment 
analysis for functional annotations of dysregulated genes. Enriched 
GO terms are statistically significant and representative in 
dysreuglated genes against the whole array. Ingenuity Pathways 
Analysis (IP A) software (Ingenuity Systems, www.ingenuity.com) 
was used to analyze IPA Canonical pathways and generate IPA 
Tox lists. Canonical pathways analysis identified the most 
significant pathways from the IPA library from the input data 
set. Tox lists were most significant toxicity lists identified from the 
IPA library [23], and were associated with key toxicological 
responses to xenobiotic insult. The p values were calculated by 
Fisher's exact test, which calculates the probability of a gene in the 
dataset of a given size due to chance alone. Since IPA software 
does not support zebrafish data, we used human orthologs of 
dysregulated genes from NimbleGen genelD file for knowledge- 
based pathway and functional analysis. Microarray data (both pre- 
and post-normalized) have been deposited in the Gene Expression 
Omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo/). 
The data is accessible through GEO Series accession number 
GSE47434 for the dose-dependent experiment, and GSE47547 
for the time-dependent experiment. 
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Results 

Mortality, growth and bioaccumulation of zebrafish 
treated with TCDD 

No significant mortality was observed in any TCDD treatment 
group after 42 d of feeding (data not shown) by Kaplan-Meier 
survival analysis with log-rank significance test (/?<0.05). 10 ppb- 
and 100 ppb TCDD caused 10% and 23% mortality in treated 
zebrafish at 42 d, respectively. Body weight and length and 
condition factor of TCDD-treated are shown in Table SI. 
Significant decreases of body weight and length were seen in 
100 ppb-treated male zebrafish at 28 d. Significant decreases of 
body length were observed in male zebrafish in the 10 and 
100 ppb TCDD-treatment groups at 28 d. Condition factor of 
female zebrafish was significantly decreased by 100 ppb at 28. 

During 42 d of feeding, accumulation of TCDD in whole body 
exhibited dose- and time-dependent increases in both female and 
male fish (Figure 3). Due to the detection sensitivity of the 
instrument, we had to pool fish samples from each group for a 
single TCDD measurement, thus there were no technical 
replicates. The amount of TCDD in female fish in the 100-ppb 
group achieved 15.49 ppb at 28 d. In male fish, the amount of 
TCDD in male fish in 100 ppb group achieved highest 
concentrations (18.04 ppb) at 28 d, but it decreased to 5.33 ppb 
at 42 d. In the 1 ppb group, there was a technical failure during 
analysis and the samples were destroyed, therefore, we lack TCDD 
assimilation data of in the 1-ppb group at 14 d. In addition, no fish 
were identified as female by physical observation at 42 d. 

TCDD-induced clinical toxicity and gross and microscopic 
lesions in zebrafish 

After 28 days of exposure to TCDD, some zebrafish (less than 
10%) in the 100 ppb group showed hemorrhaging around their 
gills. These effects became more severe at 42 days. Histopathologic 
results were characterized as three categories: epithelial lesions, 
mesenchymal lesions and edema syndrome. Male and female 
zebrafish exhibited similar degrees of lesions, which were mainly 
observed in the 10 and 100 ppb TCDD-treatment groups. Thus, 



we combined histopathologic data from both males and females, 
and the summary of lesions is shown in Figure 4. Between 28 and 
42 days, the TCDD-induced pathologic alterations exhibited time- 
and dose-dependencies in various organs of TCDD-treated 
zebrafish. 

Detailed information of epithelial lesions is shown in Table S2 
(28 d) and Table S3 (42 d). Briefly, epithelial lesions were observed 
in oropharynx (hyperplasia of esophagus and degeneration of 
pharyngeal pad), liver (depletion of glycogen), intestine (underde- 
velopment), pancreas (depletion of secretory granules in acinar 
cells), nose (underdevelopment and degeneration of neurosensory 
epithelium. Figure 5), and ovary (decreases in the fractions of 
vitellogenic oocytes compared to control fish. Figure 6). The 
evidence of TCDD-induced lesions in nasal epithelium of fish is 
novel. TCDD-induced degeneration of nasal epithelium in 
zebrafish is rarely documented. In this project, 100 ppb TCDD- 
treated zebrafish showed highly disorganized and poorly differen- 
tiated hyperplastic neurosensory epithelium compared to well 
differentiated, orderly, and uniform neurosensory epithelium in 
control fish (Figure 5). 

TCDD-induced lesions in mesenchymal tissues (Tables S4 and 
S5) included decreases in hemopoietic cells in kidney marrow 
(Figure 7), anemia in peripheral blood, foreshortened maxilla and 
diminution in abdominal adipose tissue. Edema syndromes 
(Tables S6 and S7) were observed in pericardium (Figure 8), 
posterior regions of the eyeball, and in stroma of the ovary 
(Figure 6). 

After 28 d of dietary exposure, a small percentage (less than 
20%) of 100 ppb-treated zebrafish were visually determined to be 
female by physical observation. At 42 d, none of the fish showing 
in the 100 ppb group could be determined as female by physical 
observation of the gonad. In contrast, female zebrafish in the 
control and low TCDD-treatment groups possessed matured 
ovaries, and overall, there was an equal number of males and 
females. Perinucleolar and cortical alveolar stages are two of the 
oocyte developmental stages that occur prior to the vitellogenic 
stage. In the perinucleolar stage, follicular cells are monolayer and 
oocytes are large. In cortical alveolar stages, yolk vesicles (cortical 
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Figure 3. Assimilation of TCDD in female and male zebrafish during 42-d of dietary exposure. TCDD from extraction of pooled fish (up to 
5) were analyzed by gas chromatography coupled to high-resolution mass spectrometry (GC-HRMS). The data (ppb) represent the observed TCDD 
(ng) verse wet weight offish (g). The asterisk (*) means failed analysis and no more samples to be detected in fish treated with 1 ppb TCDD at 14 d; 
the means no fish in 100 ppb group were identified as female by physical observation at 42 d. 
doi:1 0.1 371 /journal.pone.01 0091 0.gOOS 
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Figure 4. Heat-map of lesions found in dietary TCDD-treated zebrafish. Fish were collected at 28 d and 42 d after starting of dietary TCDD 
exposure. Lesions were characterized as three main categories: epithelial lesions, mesenchymal lesions and edema syndrome. Each square represents 
the percentage of fish that had the lesions. 
doi:1 0.1 371 /journal. pone.01 0091 0.g004 



alveoli) increase in size and can be observed in the follicle. Based 
on the histopathologic outcomes in Figure 6, perinucleolar and 
cortical alveolar oocytes accounted for a larger percentage in the 
ovary of the 100 ppb TCDD-treated zebrafish than in control, 
indicating that the development of the ovary in zebrafish was 
disrupted by chronic dietary TCDD exposure. Combining this 
information with sex identification by vasa expression in Table 1 
(which shows that many fish had female vasa-expression features 
even though no mature gonads were observed) further supports 
our view that development of the ovary was disrupted by exposure. 



Gene expression analysis in TCDD-treated zebrafish 

Dose-dependent experiment (0.1, 1, 10 and 100 ppb 
TCDD at 28 d). One hundred twenty-nine and 2106 genes 
were found to be differentially expressed in female and male 
zebrafish treated with TCDD at 28 d, respectively (Tables S8 and 
S9). Only 38 genes were shared in common between female and 
male zebrafish. 

Clustering analysis of gene expression patterns and top five 
enriched (p<0.05) GO terms (biological process category) for each 
cluster are shown in Figure 9. In female zebrafish, cluster 1 (44 




Figure 5. TCDD caused cystic degeneration in nasal epithelium in zebrafish after 42 d of dietary exposure. A) Well differentiated, 
orderly, uniform lamella (arrow) of zebrafish in 0 ppb group; B) cystic degeneration of nasal epithelium (arrows) was found in 100 ppb TCDD-treated 
zebrafish at 42 d (40X). 
doi:1 0.1 371/journal. pone.01 0091 0.gOOS 
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Figure 6. TCDD caused disruption of ovary development in zebrafish at 42 d. A) Ovulated vitellogenic follicles (V) in mature ovary of 
zebrafish in 0 ppb group; B) immature ovary in fish treated with 100 ppb TCDD. Ovary lacks vitellogenic follicles. Moderate interstitial edema (arrows) 
was observed; only previtellogenic oocytes are present: cortical alveolar oocytes (C) and perinucleolar oocytes (P) (10X). 
doi:1 0.1 371 /journal. pone.01 0091 0.g006 



genes) showed genes primarily down-regulated by TCDD at 
different concentrations, and these genes are mainly involved in 
lipid transport and lipoprotein assembly; cluster 2 (33 genes) shows 
genes primarily up-regulated by TCDD at different concentra- 
tions, and these are mainly involved in membrane organization 
and G-protein signaling pathways; genes in cluster 3 (52 genes) 
were up-regulated by TCDD at 0. 1 and 1 ppb but down-regulated 
by 100 ppb-TCDD, and these genes are shown to be involved in 
negative regulation of DNA replication and metabolic process. 

Dysregulated genes in TCDD-treated male zebrafish were 
divided into five clusters. Genes in cluster 1 (404 genes) were 
down-regulated in all groups except in the 10 ppb group; GO 
terms such as proteolysis, endothelial cell proliferation, steroid 
hormone receptor signaling pathway, cell proliferation and 
serotonin transport were overrepresented in these genes. Genes 
in cluster 2 (391 genes) were down-regulated in the 0.1 and 1 ppb 
groups but up-regulated in the 100 ppb group; and these are 
mainly involved in visual perception of light stimulus and sensory 
perception. Genes in cluster 3 (581 genes) were primarily up- 
regulated by 100 ppb of TCDD, and their enriched GO terms are 
related to brain development and central nervous system 
development. Cluster 4 contains 469 genes that were up-regulated 
in the 0.1, 1 and 10 ppb groups but down-regulated in the 
100 ppb group, with the enriched GO terms being myofibril 
assembly, cellular component assembly involved in morphogene- 
sis, actomyosin structure organization, striated muscle cell 
development and muscle cell development. Cluster 5 contains 
262 genes that were up-regulated at each concentration of TCDD, 
and these genes are primarily involved in cellular response to 
stimulus. These dysregulated genes were analyzed by IPA software 
and groups of genes were mapped to IPA Canonical pathways and 



Tox Lists (Table Sll). A series of IPA Tox lists was identified 
between female and male fish, such as renal necrosis/cell death, 
NRF2-mediated oxidative stress response, and cell cycle regula- 
tion; and common significant Canonical pathways were cardiac 
hypertrophy signaling, glycolysis/gluconeogenesis, p53 signaling. 
The representative top Canonical pathways and Tox lists are 
shown in Table 2. 

Time-dependent experiment (100 ppb TCDD at 7, 14 d or 
42 d). Statistically expressed genes in 100 ppb TCDD groups at 
7, 14 and 42 days are shown in Table SIO. Thirteen, 128 and 77 
genes were differentially expressed in female, immature female 
and male zebrafish, respectively, and were identified at 7 and 1 4 d 
after initiation of 100 ppb TCDD exposure; none of these were 
shared in common among groups (Figure lOA). Two hundred 
seventy-eight, 2833 and 1029 genes were found to be significantly 
dysregulated by 100 ppb TCDD at 14 d in the respective groups 
above, and only three genes were common dysregulated genes 
(Figure lOB). At 42 d, 623 and 4604 genes were significantly 
dysregulated in immature female zebrafish and male zebrafish, 
respectively; and 98 genes were common dysregulated genes 
between immature female and male zebrafish (Figure IOC). We 
mapped the differentially expressed genes to Canonical Pathways 
and Tox Lists with IPA, and statistically significant (/?<0.05) 
pathways and lists are shown in Table SI 2. Based on the Tox 
Lists, liver necrosis/ cell death, renal necrosis/ cell death and TGF- 
(3 signaling were assigned as dysregulated genes in TCDD-treated 
fish. 

Based upon physical observation and histopathologic analysis, 
none of the female fish treated with 100 ppb showed matured 
ovaries, and fraction (%) of ovaries with vitellogenic oocytes was 
near zero (Figure 6). In order to understand TCDD-induced 




Figure 7. TCDD induced hemopoietic hypoplasia in kidney at 42 d. A) Normal hemopoietic tissue (arrow) in interstitium of anterior l<idney of 
zebrafish in 0 ppb group; B) depletion of hemopoietic tissue (arrow) from interstitium of anterior kidney of zebrafish in 100 ppb group at 42 d (40X). 
doi:1 0.1 371 /journal. pone.01 0091 0.g007 
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Figure 8. TCDD caused pericardial edema in zebrafish at 42 d. A) Normal heart in zebrafish of 0 ppb group; B) pericardial edema (arrow) in 
fish fed 100 ppb TCDD at 42 d (5X). 
doi:1 0.1 371 /journal. pone.01 0091 0.gOOS 



toxicity in female fish, especially with regard to gonadal 
development, we compared the data of dysregulated genes 
between female and immature female fish using principal 
component analysis (PC A), which can evaluate differences among 
different samples based on dysregulated gene expression values. 
PCA showed a separation between mature and immature female 
fish using dysregulated genes from TCDD-treated female fish at 
each time point (Figure 11). Results from the PGA analysis 
(Figure 11) suggested that immature fish respond differently to 
TCDD than mature zebrafish. As a result, we could not 



consolidate female and immature fish into one group for gene 
expression analysis. 

We compared IPA Canonical pathways from immature fish 
with female fish; and some significant pathways were identified 
only in immature fish are related to regulation of development and 
cell death, including apoptosis signaling, growth hormone 
signaling, TGF-P signaling (Table Sll). Moreover, biomarkers 
for ovary development, such as vitellogenin 2 [vtg2), vtg4 and vtg7, 
were only found to be significantly down-regulated in immature 
fish. 



TCDD (ppb) 



10 100 



Cluster 1 
(44 genes) 



Cluster 2 
(33 genes) 



Cluster 3 
(52 genes) 




lipid transport 



lipid localization 



positive regulation of 
glycoprotein biosynthetic 
process 



plasma lipoprotein particle 
assembly 



very-low-density lipoprotein 
particle assembly 



endocytosis 



membrane invagination 



cellular membrane organization 
membrane organization 



G-protein coupled receptor 
protein signaling pathway 



negative regulation of 
biological process 



negative regulation of DNA 
replication 



ribosome assembly 



negative regulation of 
macromolecule metabolic 
process 



negative regulation of 
metabolic process 



B 



Cluster 1 
(403 genes) 



Cluster 2 
(391 genes) 



Cluster 3 
(581 genes) 



Cluster 4 
(469 genes) 



Cluster 5 
(262 genes) 



TCDD (ppb) 



10 100 




proteolysis 



endothelial cell proliferation 
steroid hormone receptor 
signaling pathway 



cell proliferation 



serotonin transport 



visual perception 



sensory perception of light 
stimulus 



sensory perception 



response to light stimulus 



cellular response to stimulus 



visual perception 



sensory perception of light 
stimulus 



sensory perception 



brain development 



central nervous system 
development 



myofibril assembly 



cellular component assembly 
involved in morphogenesis 



actomyosin structure organization 
striated muscle cell development 
muscle cell development 



phosphate transport 



cellular response to nutrient 
levels 



cellular response to 
extracellular stimulus 



cellular response to external 
stimulus 



inorganic anion transport 



0.67 



1.0 



1.5 



Figure 9. Cluster analysis (K-means) of the dysregulated genes with their representative Gene Ontology (GO) functional 
annotations from female (A) and male (B) zebrafish treated with dietary TCDD at 28 d (FDR <0.005). Data are mean fold changes for 
each treatment group. GO terms are referred to top five terms of Biological Process (BP) category. 
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Table 2. Representative significant Canonical Pathways and Tox Lists from dysregulated genes in zebrafish treated with TCDD at 
28 d. 
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Expression of exon 4 within vasa was quantified by QPCR as 
shown in Figure 12 and Figure 13. In the dose-dependent 
experiment, expression of exon 4 was significantly up-regulated in 
male fish and significantly down-regulated in females (Figure 12). 
Expression of exon 4 in males and females treated with 1 00 ppb 
TCDD were significantly different from those in control males or 
females. In the time-dependent experiment, expression of exon 4 
in males was up-regulated in both control females and immature 
female in the 100-ppb group (Figure 13); however, there was a 
sharp decrease in expression of exon 4 in 100 ppb TCDD-treated 
female fish at 28 d. 



Discussion 

Growth and accumulation of TCDD 

Exposure of zebrafish embryos (2 h post-fertilization, hpf) to 
waterborne TCDD caused decreases in total body length and 
length of bone in the jaw of lava [24]. A study by Lanham et al. 
[25] showed that body weight of zebrafish was decreased by 
waterborne exposure to TCDD at 5 day post-fertilization (dpi) or 
20-25 dpf, but decreased weight was not observed after 35 dpf, 
indicating that sensitivity of fish to TCDD exposure decreases 
during growth. In the present study, the juvenile fish were larger 
and older than those in the Lanham study; highest accumulated 
TCDD at 15.49 and 18.04 ppb in female and male fish did not 
cause significant mortality, except for decreases in body weight at 
28 d (Table SI). This indicates that life-stage and body size 
account for much of the sensitivity to TCDD in zebrafish. 





Figure 10. Venn diagram of dysregulated genes in female, immature female and male zebrafish treated with 100 ppb TCDD at 7 d 
(A), 14 d (B) and 42 d (C). 

doi:1 0.1 371 /joumal.pone.01 0091 0.gOl 0 
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Figure 1 1 . Principal component analysis (PCA) of dysregulated genes from TCDD-treated female fish at each time point. Tlie gene list 
of 3679 genes for tlie principal component analysis was generated by combining all dysregulated genes from female and immature fish at each time 
point and removing duplicates. The value for each gene within each group was computed as the ratio of average of expression of treatment relative 
to the average of control expression levels for that gene. 
doi:1 0.1 371 /journal. pone.01 0091 0.gOl 1 



Additionally, difFerences between zebrafish and other fish species 
reflect the species-specific sensitivity to TCDD exposure [26,27], 
which might be due to the specific affinities of AhR to TCDD in 
different species [28]. Decreases in body length and weight in 
TCDD-treated zebrafish may be due to loss of nutrition because of 
TCDD-induced lesions in their digestive systems (oropharynx, 
liver, intestine and pancreas) found at 28 d and 42 d. 

Previous study by King Heiden et al. [29] showed that dietary 
TCDD exposure (applied doses at 43 ppb) caused 36 ppb in 
whole-body fish at 20 d; and TCDD is highly accumulates in lipid- 
rich organs, such as ovary. In the present study, assimilations of 
TCDD in both females and males exhibited dose- and time- 



dependent increases in whole fish (Figure 3). Dietary TCDD at 
100 ppb resulted in the assimilation of 15.49 and 18.04 ppb in 
female and male fish at 28 d, respectively, without significant 
mortality. Thus, the treatments we used for exposure in this 
project can be applied to sublethal pathologic/ toxic research using 
zebrafish. 

Lesions associated with dysregulation of gene expression 
in TCDD-treated zebrafish 

Previous laboratory studies in mammals show that TCDD 
causes lesions in multiple organs, such as liver, lung [30] and 
dermis [31]. In teleost, pathologic alterations observed in zebrafish 
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Dose of TCDD (ppb) 



Figure 12. Expression of exon 4 of vasa in female and male zebrafish treated with TCDD at 28 d. The asterisk H means significant 
differences in exon 4 expression in 100 ppb group compared to control group with one-way ANOVA. 
doi:l 0.1 371/journal. pone.01 0091 0.gOl 2 
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Figure 13. Expression of exon 4 of vasa in zebrafish treated with 100 ppb TCDD within 42 d. Immature female fish were found at 7, 

14, and 42 d. And no fish in 100 ppb group were identified as female by pliysical observation at 42 d. 
doi:1 0.1 371 /journal.pone.01 0091 0.gOl 3 



were similar to those of other fish species. TCDD causes 
craniofacial malformations, disruption in jaw cartilage growth, 
yolk sac edema, altered heart rate, heat failure, and delayed 
development and death in embryos and larvae [7,29,32 35]. In 
adult fish, TCDD causes hepatic lipidosis and hypertrophy, 
glycogen depletion in liver, inhibition of fin regeneration and 
altered gonadogenesis [11,36,37]. In juvenile fish, the TCDD- 
induced lesions observed reflect similarity to exposed embryo/ 
larvae and adults, including craniofacial malformations, degener- 
ation in fin and gill, and altered pigmentation [25,38]. 

In the present study, lesions (Figure 4) identified in juvenile 
zebrafish showed certain similarities to previous research on 
juvenile yellow perch {Perca Jlavescens) treated with TCDD, such as 
decreases in hemopoietic cells in kidney marrow, depletion of 
adipose tissue, and edema syndromes [39]. Skeletal malformation 
and glycogen depletion in liver seen in the present study were also 
reported in larvae and adult fish [11,25,37]. However, altered 
pigmentation and fin necrosis were not observed, as was reported 
in 30 dpf juvenile zebrafish [25] and juvenile yellow perch [39]. 
TCDD was found to influence proliferation or differentiation 
processes of hemopoietic stem cells in bone marrow of mice after 
dietary exposure [40]. Unlike human beings, hemopoiesis occurs 
in kidney marrow of zebrafish instead of in bone marrow [41,42]. 
We found that the number of hemopoietic cells in kidney marrow 
was reduced in zebrafish treated with TCDD, and this may be 
correlated with the anemia observed in peripheral blood. 

Results from IPA showed that ahr signaling was one of 
significant Canonical pathways and metabolism of xenobiotics 
by cytochrome P450 was one of significant Tox lists (Table 2). 
TCDD is known to bind ahr and arnt to mediate gene expression, 
such as cytochrome P450, family 1, subfamily A, cypla [43]. 
TCDD-induced oxidative stress can also trigger an apoptotic 
signal and induction of cytochrome P450 enzymes is a critical 
event in the apoptotic cascade [44]. Significant enrichment 
regarding NRF2-mediated oxidative stress response, p53 signaling, 
and cell cycle regulation indicated that the lesions induced by 



TCDD in fish might correlate with TCDD-induced oxidative 
stress and apoptosis. Tox Lists related to diseases of heart, kidney 
and liver were significantly represented in both the dose- and time- 
dependent experiments, including cardiac necrosis/cell death, 
cardiac fibrosis, long-term renal injury pro-oxidative response 
panel, renal necrosis/cell death, and liver necrosis/cell death 
(Tables SI 1 and SI 2). All of these results are highly correlated with 
chronic toxic effects of TCDD we observed in heart, kidney and 
liver (Figure 4). 

TCDD-induced neuronal toxicities have mainly been reported 
in rats and mice following intraperitoneal injection of high doses. 
These include apoptosis in granule neuron cells, alterations in 
metabolism of neurotransmitters, formation of reactive oxygen 
species (ROS) in cerebellum, and dysfunction in motor, learning 
and visual abilities [45-49]. In zebrafish, embryonic exposure to 
TCDD reduces the number of neurons in the larval brain [50] and 
causes apoptosis in midbrain [51]. We reported the first evidence 
of TCDD-induced morphologic lesions in the olfactory organ of 
rainbow trout [38], in this study we found similar results from 
TCDD-treated zebrafish that exhibited highly disorganized and 
poorly differentiated hyperplastic lamellae (Figure 5). Evidence of 
TCDD-induced lesions in nasal epithelium of teleost is similar to 
results in rats which show hyperplastic and metaplastic changes 
within the nasal mucosa after chronic administration by gavage 
with TCDD [52]. 

Enriched GO terms in cluster 1, 2, 3 and 4 for treated male fish 
(Figure 9B) exhibited a correlation with nervous system develop- 
ment, such as endothelial ceU proliferation (cluster 1), brain 
development and central nervous system development (cluster 3). 
Some enriched GO terms from male fish also related to eye 
development, such as sensory perception of light stimulus and 
overall visual perception. There is a single report of TCDD- 
induced decreases in the densitiy of retinal ganglion cells in 
rainbow trout [53] but nothing in zebrafish. We did not focus on 
eye development in this study, and potential effects of TCDD on 
eye development in zebrafish need to be further explored. 
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TCDD-induced disruption of ovarian development 
associated with gene expression in zebrafish 

TCDD was found to inhibit the transition from the pre- 
vitellogenic stage to the vitellogenic stage in follicles of adult female 
zebrafish [36], and caused decreases in ovarian weight and 
ovarian necrosis [29]. Early-life waterborne exposure to TCDD 
caused impairment of reproductive capacity in both males and 
females as adults [7] . In the present study, we found that less than 
30% of the collected fish contained a mature ovary when exposed 
to 10 ppb and none in the 100 ppb group at 42 d. This made it 
difficult for us to determine their sex by physical observation 
during sampling. Histopathologic analysis showed that ovaries 
were immature in fish treated with 1 00 ppb TCDD and that many 
contained immature cortical alveolar and perinucleolar oocytes, 
stages prior to the vitellogenic stage [54] (Figure 6). In addition, 
sex determination by PCR using vasa primers also confirmed that 
the fish with abnormal gonads were female (Figure 2 and Table 1). 
Our results are thereby similar to a previous study by 
Wannemacher et al. [55] that showed significant increases in the 
number of fish with immature, previtellogenic oocytes after dietary 
exposure to TCDD. 

The vasa mRNA is expressed as two alternatively spliced 
transcripts (with or without exon 4) that are differentially expressed 
in males and females. The longer product that contains exon 4 is 
primarily expressed in females with very low expression in males 
[16]. Using our QPCR assay specifically targeting exon 4, we 
found that expression of exon 4 was down-regulated in females in 
both Dose- and Time-dependent experiments; however, expres- 
sion of exon 4 in immature females is similar to that of normal 
female fish in the control group (Figures 1 2 and 1 3), indicating that 
TCDD might only influence splicing of vasa in female fish with 
matured ovaries. Previous studies use expression of vasa mRNA as 
a marker for primordial germ cells of zebrafish at early-life stages 
[56,57], however, little is to know about the function of vasa on 
sex-determination or ovarian development during juvenile stages. 

We speculated that TCDD may inhibit the estrogen biosyn- 
thesis process by down-regulation of genes related to estrogen 
production, such as steroidogenic acute regulatory protein [star), 
cypllal, and cyplQala (aromatase), as we observed in our previous 
studies [12,36]. However, none of these genes was found to be 
dysregulated significantly by TCDD in immature fish except for 
cypl7al, which showed a four-fold up-regulation at 14 d. We also 
evaluated genes related to regulation of ovarian development, such 
as follicle-stimulating hormone receptor {fshr), luteinizing hormone 
receptor {Ihcgr), estrogen receptors [esrl, esr2a, esr2b), epidermal 
growth factor [egfj and egf receptor (egfr), and we found that only 
esrl and efgr were down-regulated at 7 d and 42 d in immature 
fish, respectively. Thus, our results cannot provide strong support 
for TCDD-induced disruption of ovarian development through 
inhibition of the estrogen biosynthesis pathway, but does suggest 
effects on receptor signaling involved in ovarian foliculogenesis. It 
is important to remember that, in this study, we performed 
microarray using RNA isolated from whole fish. In order to reveal 
transcriptional changes relevant to tissue-specific effects if may be 
necessary to use only mRNA from that tissue. 

In the Great Lakes region, TCDD levels of 0.02-0.04 ppb have 
been reported in salmonid species (e.g., lake trout) [58,59]. We 
used much higher levels of TCDD in the laboratory since we 
wished to investigate the mechanisms of dietary TCDD toxicity in 
zebrafish, which are much more resistant to the toxic effects of 
TCDD than are salmonids using other means of exposure [26- 
28,60-61]. In a recently published parallel study using rainbow 
trout [38], we found that dietary TCDD caused overt toxicity in 
juvenile rainbow trout at wide range (from 0.1 ppb to 100 ppb). 



Zebrafish in this study were more resistant to dietary TCDD 
exposure as indicated by lower observed mortality and fewer 
pathologic alterations in the lower dose groups. In conclusion, 
chronic dietary TCDD exposure from juvenile stages delayed 
maturation of ovaries in female fish and caused lesions in multiple 
organs in zebrafish, and caused dose- and time-dependent TCDD 
accumulation. Functional and pathway analysis of microarray data 
using whole fish body revealed a correlation with various lesions 
and putative mechanisms for their genesis. As for the disruption of 
ovarian development and lesions in nasal epithelium induced by 
TCDD, tissue-specific microarrays should be performed to explore 
the mechanisms. 
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